Superconducting Dome from Holography 
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We find a regime in which a strongly coupled striped superconductor features a superconducting 
dome. This regime is signified by i) a modulating chemical potential that averages to zero, and ii) 
a superconducting order parameter that has a scaling dimension 3/2 < A < 3. We also find that in 
this regime, the order parameter exhibits a mild dependence on the modulation wavelength of the 
stripe. 
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Quantum many-body physics lacks a general mathe- 
matical frameworks that can deal with strongly coupled 
fermions at finite density. In an effort to remedy this 
problem, recently, there has been a flurry of activities 
applying gauge/gravity duality to study strongly cou- 
pled condensed matter systems (for a review, see, e.g., 
PQ). A large portion of these activities revolves around 
understanding high temperature (high T c ) superconduc- 
tors, such as cuprates and iron pnictides. One of the 
many challenges in adopting this approach is to repro- 
duce the phase diagrams of the high T c superconductors. 
The problem is that it is not clear how one should intro- 
duce the effect of doping into the holographic model. 

Being strongly coupled, the normal states of high T c 
superconductors are highly correlated and thus, exhibit 
other low temperature orders which interact with su- 
perconductivity in an intricate way. One of the promi- 
nent orders is the unidirectional charge density wave or 
"stripe" order [2H5]. However, the complication that 
comes from the existence of these coexisting orders might 
also be the key toward reconstructing the phase diagrams 
of high temperature superconductors. In particular, a 
very important hint comes from the observation that in 
cuprates, in particular Bi2_yPbj,Sr2_ z La z Cu06 +a; , near 
the superconducting regime the modulation wavenum- 
ber is a monotonically increasing function of doping [5] . 
Because of this, the qualitative behavior of the critical 
temperature viewed as a function of either doping or the 
modulation wavenumber should be similar. Therefore, as 
a step toward realizing a realistic holographic model of 
cuprates, it is crucial to show that indeed there exists a 
regime in the parameter space in which the holographic 
striped superconductor exhibits a superconducting dome. 

To this end, following [6j, we consider an Einstein- 
Maxwell-scalar system on a 3 + 1-dimensional spacetime 
with negative cosmological constant A = — 3/L 2 as a 
holographic model of a 2 + 1-dimensional strongly cou- 
pled striped superconductor. The bulk scalar field is dual 
to the superconducting order parameter, while the bulk 
U(l) gauge field is dual to the four-current in the strongly 
coupled system. For simplicity, we shall adopt units in 



which L = 1, IGirG — 1. We introduce the stripe order 
phenomenologically by introducing a modulating chemi- 
cal potential ia[x) = n(l — 5) + fi5 cos Qx, which is identi- 
fied as the boundary value of the bulk electrostatic poten- 
tial. We are interested in a black hole solution with this 
boundary condition and the superconducting transition 
temperature will then be given by the critical tempera- 
ture below which the black hole forms scalar hair. 

It has been shown that if the average value of the chem- 
ical potential does not vanish [6 ^ 1), then the critical 
temperature will have a long power law tail at large Q. 
This is true both at the mean field leveQ [7] and when 
fluctuations are turned orj^] [6! . A similar result is also ob- 
tained when 5=1, and the scaling dimension of the order 
parameter satisfies 1/2 < A < 3/2. This is an unwanted 
feature since to obtain a superconducting dome, we need 
to have a critical modulation wavenumber > 0, above 
which the critical temperature vanishes. At the mean 
field level, it was shown in Ref. [7 that the 5 = 1 and 
3/2 < A < 3 case results in such a feature. 

Another hint we obtained from our previous work is 
that in the 5^1, and the 5 = 1, 1/2 < A < 3/2 cases, 
fluctuations cause a significant drop in the critical tem- 
perature at small Q [B]. Combining the above, we expect 
that when the chemical potential averages to zero (5 = 1) 
and the scaling dimension of the order parameter satisfies 
3/2 < A < 3, the superconducting regime is capped at 
larger values of Q, while when the fluctuations are turned 
on, it is capped at the smaller values of Q, thus resulting 
in a superconducting dome. We will show that indeed 
this expectation is correct. 

The equations of motion consist of the Einstein equa- 
tions, R a b — 9abR/2 — igab = Tab/2, the Maxwell equa- 
tions, d b (xf I gF ab )/ x f^g~ = J a , and the Klein-Gordon 
equation for the scalar field, —D a (y/—gg ab Db<fi) j \J— g + 



1 This is equivalent to the probe limit in the gravity picture, where 
the backreaction of the U(l) field and the scalar field on the 
spacetime geometry is neglected. 

2 Equivalently, when backreaction is included. 
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m 2 (f) = 0. Here, T a b is the stress-energy tensor, the 
U(l) current is J a = —i[<fi*D a (f) — c.c.], the covariant 
derivative is D a — d a — iqA a , the U(l) field strength is 
F a b = d a Ab — dbA a and a,b £ {t, z, x, y}. Furthermore, 
q is the charge of the scalar field, which is related to the 
central charge of the conformal field theory describing 
the critical point and m is the mass, which is related to 
the scaling dimension A by m 2 = A(A — 3). 

We solve the equations of motion by doing an expan- 
sion in 1/q 2 , where taking q — ¥ oo while keeping /i fixed 
corresponds to the probe limit. In the normal state, 
where the order parameter vanishes and therefore, = 
, the solution for the metric for arbitrary S is given in [6] . 
Up to 0(l/q 2 ), the solution for 5 — 1 is then given by 
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where the boundary is at z = and the horizon is at 

a^{z)+a^ 2 \z) cos2Qx, 
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Prime ' denotes a derivative with respect to z. The tem- 
perature is given by 
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The electrostatic potential is given by A t 
[Aq + A\/q 2 ) cosQx, where Aq and A\ satisfy 
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As we lower the temperature, = becomes an un- 
stable solution to the scalar equation of motion 
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To solve this, we expand in a Fourier series 
(O a )z a 
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where (Oa) is the superconducting order parameter. We 
evaluate the critical temperature numerically and a typ- 
ical result can be seen in Fig. [TJ 
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FIG. 1: (Color online) The critical temperature as a func- 
tion of modulation wavenumber. Here, A = 3 and the 
black, red and green lines correspond to q 2 — oo, 10 and 
5, respectively. 



Turning on fluctuations is equivalent to decreasing the 
value of q 2 and this results in the critical temperature 
decreasing at small Q. The discontinuous behavior ex- 
hibited at Q — > is an artifact of the first-order con- 
tribution to the perturbative expansion of the solutions 
of the Einstein-Maxwell equations. We expect that as 
we include higher orders in perturbation, the graph will 
become smoother. 

In the regime of validity of the 1 j q 2 expansion, we are 
not able to crank up the fluctuations to reach vanishing 
critical temperature at Q = 0. However, as the exact 
solution at Q = is known, we expect this to happen at 
q 2 = 3/4 + A(A-3)/2 0. 

Unlike the cases considered in Ref. [Hj, here we see 
that there is a critical modulation wavenumber , above 
which the critical temperature vanishes. The value Q* 
can be estimated analytically by considering the large 
modulation wavenumber regime where Q > T c . In this 
regime, we can see from the numerics that the effects 
of backreaction are suppressed. Furthermore, the higher 
modes are suppressed. The equation of motion is then 
given by 



d 2 F (o) + 2(A-l) d _ pW 



+ Xe- S F^ = 0, (8) 



where z = zr+/(2Q) and A = q 2 fi 2 /(8Q 2 ,). We would 
like to evaluate this in the interval z 6 [0, oo), with the 
boundary condition — > as z —> oo. This differential 
equation can be solved exactly for A = 2, but for other 
values 3/2 < A < 3, we can estimate A by iteration. To 
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FIG. 2: The end point of the superconducting dome as 
a function of the scaling dimension of the order 
parameter. 



do so, let us rewrite [8] as an integral equation 
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Solving this iteratively, we find that at zeroth order 



~qil) ~ 8(2A-3)' 
while a hrst order correction results in 
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Here, ^ denotes the digamma function. Despite appear- 
ances, as we can see in Fig. [2j this is a smooth function 
of A in the interval (3/2,3) (the two expressions above 
have almost indistinguishable graphs). Furthermore, at 
A 4 3", we have ^* ~ 0.158, in good agreement with 
our numerical results. 

By comparing the information concerning the end 
points of the dome, Q — and Q — Q*, with experi- 
ment, we can then in principle extract the value of A. 
This is because the Q = point gives us the information 



on q and upon substituting it into Eq. 11 we obtain A. 

Another method to extract the realistic value of A is 
by calculating the anisotropy of the optical conductivity 
akin to the calculation done in [5] (see also [TO]) and com- 
paring its scaling behavior with the observed behavior in 
cuprates |llj . Work in this direction is in progress. 
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FIG. 3: (Color online) The ratio of the gap with respect 
to temperature as a function of the wavenumber at fixed 
T/T c . Here, A = 3 and the blue, purple and yellow lines 
(bottom to top) correspond to T/T c = 0.975, 0.95 and 
0.9, respectively. For clarity, we have also shown the 
critical temperature, plotted in dashed line. 



Now that we have understood the behavior of the sys- 
tem at the critical temperature, we can go to a tem- 
perature T below T c and study the behavior of the order 
parameter. Limiting ourselves to the case where the bulk 
scalar field can be considered as a small perturbation to 
the solution at critical temperature, we obtain the gap 
at fixed T/T c as plotted in Fig. [3j This corresponds to 
the temperature below but near T c . 

We see that deep under the superconducting dome, the 
gap shows a mild dependence on Q. Interestingly, mea- 
surement of the gap at fixed temperature at low temper- 
ature also showed mild dependence on Q [T^]. We are 
hopeful that gap measurement near T c from STM will be 
available to be compared with our prediction in Fig. [3] in 
the near future. 

The ratio (Oa) 1 ^ A /T deep under the superconducting 
dome as a function of A is plotted in Fig. [4] The value 
is about 8-10, which is the same order of magnitude with 
the measured value of 29 [5] . We would like to note that 
the measurement is done at T/T c ss 0.2, which is beyond 
the regime of validity of our approximation. However, 
considering that in some holographic models, the ratio of 
gap over temperature can go as high as O(10 3 ) (see for 
example [13 ) , our result is encouraging. 

The result at the far end of the dome Q ~ Q* can be 
understood analytically using the perturbative method 
of Ref. 0. We have 
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FIG. 4: (Color Online) The ratio of the gap to 
temperature deep under the dome, as a function of the 
scaling dimension of the order parameter. Here, the 
blue, purple and yellow lines (bottom to top) 
correspond to T/T c = 0.975, 0.95 and 0.9, respectively. 



where 



a(z) 



a(z) 



I z 



dz ' F (0) ^£ p(0) 
2-2A r h r > 



o r 



dz 1 

/2-2A 



Here. 

F (0) 

and 



r , A A 2A t 

tFx[ j>V~3~'' z 



(13) 



(14) 



_3-2A v 



3 - A 3 - A 2(3 -A) 



-2Qz/r. 



V 2A-1 



1 



(0)' 



(15) 



2 ^ (2A-1)Q r^- 1 

We note that by setting the order parameter to zero in 
Eq, ~ 
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we can obtain the value for r\J (<? 2 /i 2 ) from Eq. 
12 where r +c is the value of r + at the critical temper- 
ature. Then, by dividing Eq. 12 with r 2 +c /(q 2 fj 2 ), after 
some algebra we obtain 
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We see that as Q — > Q* that even though the gap van- 
ishes, the ratio of the gap to critical temperature diverges 
as T c 3 ^ 2A \ It is interesting that we found a regime in 
which the value of (O a ) 1/a /T c gets very large at low 
temperature. In comparison, the values for homogeneous 



BCS superconductor and homogeneous holographic su- 
perconductor are 3.54 and 8, respectively [13]. 

Summary and Discussions - In this article, we have 
shown that when the homogeneous part of the chemi- 
cal potential vanishes and when the scaling dimension of 
the order parameter is larger than 3/2, the holographic 
striped superconductor exhibits a superconducting dome 
with vanishing or nearly vanishing critical temperature 
at the end points. Fitting the features of the end points 
to experimental data will give us the central charge of 
the conformal field theory describing the critical point 
and the scaling dimension of the superconducting order 
parameter. 

We have also shown that deep under the supercon- 
ducting dome, the order parameter exhibits a mild de- 
pendence on the modulation wavenumber. This is per- 
haps somewhat unexpected and shows that the lack of 
Q-dependent features in the superconducting gap data 
should not be taken as an indication that the stripe or- 
dering is not related to the mechanism that gives rise to 
superconductivity. 

Toward the end of the superconducting dome, the gap 
exhibits qualitatively different behavior. In particular, as 
one approaches Q = Q*, the gap decreases to zero while 
its ratio with respect to critical temperature increases. 

Lastly, we would like to remind the reader that in this 
holographic model, the stripe order is introduced phe- 
nomenologically in the form of a modulated chemical po- 
tential. In order to understand the phase diagram of the 
holographic striped superconductor more accurately, we 
ultimately would like to have both the stripe order and 
superconducting order emerging dynamicalljj^J Work in 
this direction is in progress. 
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